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Abstract 

Strutinsky shell corrections for the Cesium-coated fullerenes were investigated. The sin- 
gle particle levels of electrons are obtained using the spherical mean-field potential of a 
shifted Wood-Saxon type. The parameters of the potential are adjusted to reproduce the 
experimental ionization energies of the Cs(iV) clusters and the magic numbers observed 
in their photo-ionization spectra of the C6oCs(iV) aggregates. 

PACS. 36.40.Cg Electronic and magnetic properties of clusters 71.24.+q Electronic struc- 
ture of clusters and nanoparticles 



1 Introduction 

The world wide interest on metallic clusters started in the early 80's when a group of physicists 
from Berkeley performed the first experiments with tiny alkali clusters [|l], |2|. Since that time an 
enormous amount of work both on the experimental and theoretical side has been accomplished. 
Metallic clusters are the quantum many body systems built from several up to a few thousands 
atoms and their radii are of the order of 10 A. It is well known that the shell structure of 
electronic orbitals in the metallic cluster can be well reproduced not only by selfconsistent 
calculations but also within some average potential approach based on the jellium model (see 
e.g. 0). In the jellium model one assumes, that electrons are moving in the mean field potential 
which is produced by the background of positive ions and by other free electrons. 

In 1996 an experimental group from Stuttgart || succeeded in producing and the measuring 
optical response of Cesium coated fullerenes Ceo covered by a layer of iV Cs-atoms with N < 500. 
The electronic shell structure of these C6oCs(iV) agglomerates was studied in Ref. using 
the local approximation to the density-functional theory. The measured abundances of metal- 
coated fullerenes C6oCs(iV) were interpreted there in terms of the magic numbers corresponding 
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to the electronic shell closures. Recently, in Ref. the shell structure of thin spherical metallic 
layers was discussed using the Strutinsky shell correction method |7[] and an infinite square 
well approximation to the mean field potential for electrons. In spite of this rather rough 
approximation to the average mean field potential it was shown there that for C6oCs(iV) the 
minima of the shell correction as function of the atom number N in the cluster were found not 
far from the experimental magic numbers given in Ref. [||]. The agreement was better at higher 
mass numbers, where the jellium model is more reliable. 

The aim of the present investigation is to use a similar approach as in Ref. || but with a 
more realistic approximation to the mean-field potential. We are going to keep the spherical 
symmetry of the potential because the fullerene Ceo core is a spherical ball, what additionally 
stabilizes the system against deformation. We show in the following, that a shifted Saxon- 
Woods potential reproduces reasonably well the experimental magic numbers and the ionization 
energies. It will also be shown how the magic numbers of the Cesium layer change with the 
core radius. This is not only an academic study since one can easily imagine an experiment in 
which one covers a spherical insulator with a thin metal layer. 



2 Theoretical model 

The valence electrons have a decisive influence on the abundance of metallic clusters. The 
electrons are moving in a mean field potential which is produced by the jellium, i.e. the 
positively charged ions and the rest of the electrons. 

It was shown (see discussion in Ref. 0) that the Saxon- Woods type potential 

V phen (r) = , (1) 

approximates quite well the mean field which is felt by free (valence) electrons. The radius Rq 
in (JTJ) has to be chosen as 

R = roiV 1 / 3 , (2) 

with 

a 

r = r 00 (l + j^) , (3) 

where r 00 is the Wigner-Seitz cell radius, N the number of atoms in the cluster and a a 
constant originating from the so-called "spilled-out" effect [1], [|. The potential depth V and 
the surface thickness a are usually adjusted to the experimental values of ionization energies 
and to some spectroscopic data (e.g. magic numbers) of clusters in a broad range of N. Solving 
the Schrodinger equation 

(~2^r y2 + Vphen ) ® u = e ^ u ' ^ 

one can obtain the ionization energy by looking for the energy of the last occupied level (Fermi 
energy) . One has to keep in mind that due to the spin of the electron one has to put two electrons 
on each single-particle level e v . In the spherical case the levels are strongly degenerated 

^u = ^n£mm s , while C v = e ni . 

The magic numbers can be roughly identified with the numbers of particles below energy gaps 
in the single-particle spectrum e u or more accurately as minima in the considered as function 
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of the particle number N Strutinsky shell correction 5E sheU (N) |7|, |8| 

6E ahea (N) = 2 J £e u -E(N) , (5) 

occ 

where E is energy of electrons in the cluster with smeared-out shell structure. We have used 
here, as usually done in the Strutinsky procedure, a Gauss function with the six order correc- 
tional polynomial to washed out the shell effects. The width (7) of the Gaussian was fixed from 
the plateau condition of 

00 

E(N) = 2Y,e»n v (l) , (6) 

where n u is the average occupation number of the level e„. 

We are going to discuss the shell structure of the electron orbitals in thin Cesium layers 
built on a fullerene core or another spherical ball made from insulator. A shifted Saxon- Woods 
potential seems to be most appropriate in this case. Such a potential, shown in Fig. 1, is 
described by the following equation 

V(r) = V ( ] —\ , (7) 

\l + e a 2 l + e a i / 

where R2 is the radius of the inner core and Ri depends on the number of atoms in the cluster 

R 1 = (r 3 N + R 3 2 ) 1 / 3 . (8) 

The parameters a\ and a 2 describe the diffuseness of the outer and inner potential wall at R\ 
and i? 2 respectively. The parameter a x should be essentially equal to the parameter a from 
Eq. (0) for the compacted cluster. The potential which feels the electron in the region between 
the core and the metallic layer of the Cesium atoms is described by the parameter 02 which 
will be adjusted to reproduce the observed magic numbers in the CeoCs(iV) agglomerates as 
explained in the next section. 



3 Results 

The single-particle levels of electrons were evaluated using the Woods-Saxon potential (|l|) for 
the compact spherical clusters build from a few to 1000 atoms. The magnitude of the ionization 
energy dependens on the depth V$ of the potential, its radius Ro, and surface thickness a: 

E io = e N / 2 {Vo,Ro,a) . (9) 

In order to fix the parameter set of the potential (|I|) we have performed a preliminary 
computation for Sodium Na(iV) and Potassium K(N) clusters which are the best known clusters 
of metals with one valence electron and for which the ionization energies were measured in 
Refs. ||, in a relatively broad range of N. 

We have found that using the radius shift constant a = 0.8 and the depth of potential 
Vo = 6.0 eV for Sodium and 4.6 eV for Potassium clusters, and the surface thickness a=0.3 A 
for both type of clusters one can reproduce quite accurately the experimental systematics of 
ionization energies as well as their average phenomenological dependence on [|T(| : 

E - en (N) = l^TT 3 -EUoo), (10) 
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where E io (oo) is the ionization energy equal to 2.7 eV for Sodium and 2.4 eV for Potassium 
metals. 

The comparison of the theoretical and the experimental ionization energies and their average 
dependences on iV is presented in Fig. 2 for the Sodium clusters and in Fig. 3 for the Potassium 
clusters. The experimental magic numbers for Sodium and Potassium as well as the other 
alkaline clusters are the same and equal 

N exp = 2, 8, 18, 20, 34, 40, 58, 92, 138, 198, 264, 344, 442 . (11) 

The theoretical magic numbers closest to the experimental ones 

Nfh = 2, 8, 18, 20, 34, 40, 58, 92, 138, 186, 198, 254, 338, 440 , (12) 

which are obtained with the diffuseness parameter of the Saxon- Woods potential equal to a=0.3 

A. 

The Cesium with the ionization energy £'j (oo)=3.89 eV belongs also to the group of the 
alkaline metals, so the Cs clusters should have similar properties as the Na and K clusters. It is 
the main reason why we have adopted the same value of the radius shift parameter (|J) a = 0.8 
performing calculations for the Cs aggregates. The theoretical systematics of the ionization 
energies for the Cs(iV) clusters obtained with the depth of the Saxon- Woods potential equal 
to Vq = 5.6 eV and with the surface thickness a=0.3 A is very close to its phenomenologic 
estimate ]10| as it can be seen in Fig. 4. If the surface thickness a is larger than 0.5 A the 



asymptotic behavior of the ionization energy differs significantly from the estimate made in 
Ref. ||1Q|| . In Fig 5 the Strutinsky shell correction for the compact spherical Cesium clusters is 
plotted. The minima of the shell correction (|5|) correspond to the predicted magic numbers. It 
is seen that the theoretical magic numbers are exactly the same as the experimental ones for 
N < 138 and very close to the measured ones for heavier clusters. 

Varying the inner and outer surface widths of the shifted Saxon- Woods potential ([TJ) and 
keeping the rest of the parameters constant we have tried to find a combination of the a% 
and a,2 widths leading to the magic numbers experimentally found in Ref. HQ for the CgoCs(iV) 
aggregates. The experimental and theoretical magic numbers from Ref. H obtained in the local 
density-functional theory are shown in the first two columns of Table 1. The magic numbers 
obtained in Ref. |1 for the square well potential are printed in the third column. In Columns 
4-9 we have displayed the magic numbers obtained with the potential (|7|) for different sets of the 
surface thicknesses. We begin with the outer surface thickness ai significantly larger than the 
inner one (02) in Column 4 and we end with an opposite choice of the thicknesses in Column 9. 
It is seen that, in spite of the dramatic change of the parameters, the predicted magic numbers 
for the lighter clusters do not differ significantly. More pronounced differences are found for the 
middle mass and heavy Cs clusters. We follow the suggestion made in Ref. 5 when evaluating 
the positions of the magic numbers in C6oCs(iV) and assume that 6 electrons are accumulated 
in the fullerene core. 

There is no reason why the parameter a\ in Eq. (|7|) for CgoCs(A^) should not be equal to 
the parameter a in Eq. ([!]) for the compact cluster. So, in the following, we have restricted 
our model and we keep a\ = 0.3 A and we change only the inner surface thickness. The magic 
numbers closest to the experimental data obtained with the above restriction are shown in the 
last column of Table 1. An example of evaluated shell corrections for C6oCs(iV) clusters is 
plotted in Fig. 6. The minima of 5E s h e ii correspond to the magic numbers predicted by our 
model when one take into account the electrons being inside the Cqo core. One can observe 
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some differences in the positions of the minima in comparison with Fig. 5. This is due to the 
effect of the inner potential wall originating from the fullerene core. The amplitude of the shell 
correction reaches 1 eV for the heavier C6oCs(iV) clusters. 

The electron single particle levels in the C 60 Cs(iV) aggregates are plotted in Fig. 7 as function 
of the thickness of the Cs layer which is equal to the difference between outer radius R\ and 
the radius of the fullerene core R<i- The positions of the Fermi level for different cluster sizes 
are marked by the black points. It is seen in Fig. 7 that for the large thickness of the Cesium 
layer the effect of the repulsive core in the single particle potential (^) becomes negligible while 
for the smaller thicknesses the level corresponding to the low angular momentum are pushed 
up. This figure is very similar to the single particle plots presented in Ref. || for the metallic 
bubbles. 

4 Summary 

The phenomenological Woods-Saxon potential with an repulsive inner core describes relatively 
well the values of the magic numbers observed in the abundance of the CeoCs(N) clusters. We 
have shown how strongly the structure of the single-particle electron levels changes with the 
thickness of the metal layer covering the spherical core. 

The results of our investigation could be also important for practical or even technological 
applications, because we indicate how one can tune in an almost continuous way the eigenen- 
ergies of electrons. In other words we give a prescription how to change the optical response of 
thin metallic layers covering arbitrary spherical cores made e.g. from insulator or other material 
with significantly different properties, by changing the thickness of the layer or the core size. 

References 

[1] W. D. Knight, K. Clemenger, W. A. de Heer, W. A. Saunders, M. Y. Chou, M. L. Cohen: 
Phys. Rev. Lett. 52, 2141 (1984) 

[2] W. D. Knight, W. A. de Heer, K. Clemenger, W. A. Saunders: Solid State Commun. 53, 
44 (1985) 

[3] M. Brack: Rev. Mod. Phys. 65, 677 (1993) 

[4] S. Frank, N. Malinowski, F. Tast, M. Heinebrodt, I.M.L. Billas, T.P. Martin: Z. Phys. 
D40, 250 (1997) 

[5] M. Springborg, S. Satpathy, N. Malinowski, T.P. Martin, U. Zimmermann: Phys. Rev. 
Lett. 77, 1127 (1996) 

[6] K. Pomorski, K. Dietrich: Eur. Phys J. D4, 353 (1998) 

[7] V.M. Strutinsky: Yad. Fiz. 36, 614 (1966) 

[8] S. G. Nilsson, C. F. Tsang, A. Sobiczewski, Z. Szymahski, S. Wycech, G. Gustafson, I. L. 
Lamm, P. Mdller, B. Nilsson: Nucl. Phys. A131, 1 (1969) 

[9] W. A. de Heer, Rev. Mod. Phys. 65, 611 (1993) and references therein. 



5 



[10] A. Tamura: Eur. Phys. J. D9, 249 (1996) 



6 



Table 1. Comparison between the magic numbers for several theoretical and experimental 
results. The single particle potential thickness parameters ai and a 2 are in A units. 
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Figure captions: 

1. The shape of the spherical average potential for electrons in the C6oCs(iV) cluster. 

2. The ionization energy of Sodium clusters Na(iV). 

3. The ionization energy of Potassium clusters K(N). 

4. The ionization energy of Cesium clusters Cs(iV). 

5. The shell-correction energy and the magic numbers for the Cs(iV) cluster evaluated with 
the surface thickness a=0.3 A. 

6. Shell-correction energy and the magic numbers for the C6oCs(iV) cluster evaluated with 
the surface thickness ai=0.3 A and a 2 =0.1 A. 

7. The single-particle levels in the Cesium layer covering the spherical core of radius R 2 - 
The depth of the single particle potential is equal to 5.6 eV, and the surface thicknesses 
ai=0.3 A and a2=0.1 A. The black points show the position of the Fermi level. 
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Figure 1: The shape of the spherical average potential for electrons in the C 60 Cs(./V) cluster. 
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Figure 2: The ionization energy of Sodium clusters Na(iV). 
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Figure 3: The ionization energy of Potassium clusters K(N). 
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Figure 4: The ionization energy of Cesium clusters Cs(N). 
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Figure 6: Shell-correction energy and the magic numbers for the C 60 Cs(iV) cluster evaluated 
with the surface thickness ai=0.3 A and a 2 =0.1 A. 
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Figure 7: The single-particle levels in the Cesium layer covering the spherical core of radius 
R 2 - The depth of the single particle potential is equal to 5.6 eV, and the surface thicknesses 
ai=0.3 A and a 2 =0.1 A. The black points show the position of the Fermi level. 
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